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Quenched QCD at finite temperature with overlap Fermions * 
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We study quenched QCD just above the phase transition temperature using overlap Fermions. Exact zero 
modes of the overlap operator are localized. Chiral symmetry is restored, as indicated by the behavior of the 
chiral condensate after subtracting the effects of zero modes. The vector and pseudo-scalar screening masses are 
close to ideal gas values. 



In QCD for T > T c screening correlators in the 
vector, and axial- vector channels are saturated by 
nearly non-interacting quark anti-quark pairs in 
the medium On the other hand, the scalar 

and pseudo-scalar screening masses show more 
complicated behavior — strong deviations from 
the ideal Fermi gas, and a strong temperature de- 
pendence. This puzzling behavior has been seen 
in quenched || and dynamical simulations with 
two Q and four flavors (l]j||| of staggered quarks, 
as well as with Wilson quarks Q close to T c . The 
new technique we bring to bear on this problem is 
the use of overlap Fermions [Q . It has the advan- 
tage of preserving chiral symmetry on the lattice 
for any number of massless flavors of quarks || , 
and hence having the correct number of pions. 

The overlap Dirac operator (D) can be defined 
JqJ in terms of the negative Wilson-Dirac operator 
(D w ) by the relation 



D = l-D w (DlD w ) 
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(1) 



The computation of D^ 1 needs a nested series 
of two matrix inversions for its evaluation; each 
step in the numerical inversion of D involves the 
inversion of D^-D^. Details of the computation of 
the inverse square root of the matrix can be found 
in The study was performed on quenched 

QCD configurations at temperatures of T/T c = 
1.25, 1.5 and 2 on 4 x 8 3 and 4 x 12 3 lattices. The 



corresponding couplings are respectively /3 = 
5.8941 and 6.0625. 

A massive overlap operator is defined by 
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D(ma) = ma+(l-ma/2)D = G~ 1 (ma), (2) 

where m is the bare quark mass, a the lattice 
spacing, and D is defined in (JIJ). We computed 
the quark propagator, G, on 12 point sources 
(3 colors and 4 spins) for 10 quark masses from 
?7ia=0.001 to 0.5 using a multimass inversion of 
D^D. The Wilson mass term in D Wl which is 
an irrelevant regulator, was set to 1.8. The tol- 
erance was e = 10~ 6 in the inner CG and 10 -4 
in the outer CG. This meant that the Ginsparg- 
Wilson relation |ll[] was satisfied to an accuracy 
of 10~ 9 , and a chiral Ward identity to an accuracy 
of 10~ 5 . 

Seperately, a rough computation of the eigen- 
values of D^D, /i 2 , was made on each configura- 
tion with a Lanczos method in each chiral sector. 
Whenever /i 2 ~ O(10~ 5 ) was obtained, the few 
lowest eigenvalues and eigenvectors were refined 
to a precision of about 10~ 8 by a Ritz functional 
minimization. As shown in Figure [l], for most 
configurations we found that the lowest /x 2 was 
well away from zero. However, for some configu- 
rations we found zero and near-zero modes with 
/i 2 < 10~ 4 clearly seperated by a gap from the 
non-zero modes with 0.1 < /x 2 . The non-zero 
modes clearly came in degenerate pairs of oppo- 
site chiralities. The zero modes were all less than 
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10 -7 and of definite chirality. There seems to be 
a spectral gap between these and the near zero 
modes (which were seen only at T = 1.25T C ). It 
remains to be seen whether the gap scales with 
lattice volume. 
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Figure 1. Scatter plot of the eigenvalue, fi 2 , 
and localisation, a for the smallest eigenvalues 
on 4 x 12 3 lattices at T = 1.25T C . The pluses and 
crosses denote chiral positive and negative sectors 
respectively. 



We constructed a gauge-invariant measure of 
localization, a |l2[ , which varies from unity for 
an eigenmode localized at just one site to 1/V 
for an eigenmode spread uniformly on a lattice of 
volume V. As shown in Figure |l|, zero and near- 
zero modes tend to be more localized than the 
non-zero modes. 

Our measurement of (tpip) comes from the diag- 
onal part of G(ma) on the 12 sources used for each 
configuration. Writing G in terms of the eigenvec- 
tors of D'D with eigenvalue [i 2 and chirality 
a, the contribution of a zero mode can be eas- 
ily read off from the equation above and is seen 
to be proportional to <f><I>/ma. Since the eigen- 
vector corresponding to the zero mode is local- 
ized, its contribution to the condensate depends 
strongly on the spatial position of the source vec- 



tor. The remaining modes are delocalized and 
closely spaced; so the overlap of the eigenvector 
on the source is averaged out. After subtract- 
ing the zero mode contribution, (iptp) is strikingly 
identical to that seen in the sample without zero 
modes, at all the couplings and lattice sizes stud- 
ied. We found that (ipip) varies linearly with ma 
and goes to zero as ma — > 0, with a value of 
(ipTp)/ma which is independent of the lattice vol- 
ume. This is how chiral symmetry restoration 
manifests itself in quenched QCD. 
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Figure 2. a 3 (tpi/j)/ma as a function of ma at 
T = 1.25T C (circle), 1.5T C (diamond) and 2T C 
(square) on 4 x 12 3 lattices. Zero and near- 
zero mode contributions are subtracted. Data for 
1.5T C and 2T C are displaced in ma for visibility. 



The following identities hold for overlap 
Fermions in the chirally symmetric phase as 
ma — > 0, 

C s (z) = -Cps(z) and C v (z) = C AV (z). (3) 

Here C is the screening correlation function in the 
spatial z direction of an operator summed in the 
other three directions. The subscripts PS refer 
to a pseudo-scalar operator, S to a scalar, V to a 
vector and AV to an axial-vector. We find that 
the V and AV correlators indeed agree at all z and 
at all temperatures we studied. The S and PS 
obey this relation after the zero modes have been 
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subtracted. In addition, Cy is described well by 
an ideal gas of overlap quarks on the same lattice, 
while Cps gives a screening mass within 10% of 
the ideal has result (see Figure || for an example). 




Figure 3. The screening correlators on 4.12 3 lat- 
tices at T = 1.5T C for ma = 0.001, in analy- 
ses without zero modes. The V/AV correlators 
(boxes) agree very well with the ideal gas compu- 
tation (solid line), and the S/PS correlators (cir- 
cles) are also similar. 



In conclusion, working with chiral (overlap) 
Fermions, we have found several new results and a 
consistent picture of the high temperature phase 
of quenched QCD. The quenched thermal en- 
semble contains gauge fields which give rise to 
Fcrmion zero modes of definite chirality. When 
the effect of these modes is subtracted, (i/ji/j) van- 
ishes in the zero quark mass limit, showing that 
chiral symmetry is restored. Simultaneously, par- 
ity doubling is seen in the spectrum of screening 
masses, which are close to those expected in an 
ideal Fermi gas, even for the S/PS sector. Since 
some of these results are not obtained with stag- 
gered quarks, it is an interesting question whether 
the two flavor QCD phase transition is properly 
described by such a representation of quarks. 

Some interesting problems remain to be solved. 
At T < 1.25 T c , there are near-zero modes. 
It cannot be ruled out that these modes shift 



the quenched chiral symmetry restoration point 
away from T c . However, this question is cru- 
cially related to the evolution of near-zero modes 
with lattice volume and spacing. Hence the na- 
ture of these complications, and the question of 
whether they are quenched artifacts or remain 
in full QCD, will only become clear with further 
studies which are underway. 
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